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Frustrated magnets such as spin ice exhibit Coulomb phases, where correlations have power-law 
forms at long distances. Applied perturbations can cause ordering transitions which cannot be 
described by the usual Landau paradigm, and are instead naturally viewed as Higgs transitions of 
an emergent gauge theory. Starting from a classical statistical model of spin ice, it is shown that 
a variety of possible phases and transitions can be described by this approach. Certain cases are 
identified where continuous transitions are argued to be likely; the predicted critical behavior may 
be tested in experiments or numerical simulations. 
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I. INTRODUCTION 

The frustrated magnetic compounds known as spin 
icejliJ^ including Ho2Ti207 and 073X1207, are well de- 
scribed by a classical Ising model on the pyrochlore lat- 
tice. As originally noted by Anderson,!^"^ this model has 
an extensively degenerate low-energy manifold closely 
related to that of water icej^ experiment^ indeed ob- 
serve the characteristic low-temperature residual entropy. 
Neutron scattering reveals long-ranged dipolar correla- 
tions in these materials, which are in close agreement 
with numerical simulations. - 

Such power-law correlations are a signature of the 
"Coulomb phase"!^ common in frustrated systems, where 
the low-energy manifold is highly degenerate but strongly 
constrained. Models such as spin ice, where degrees 
of freedom on the links of a lattice are constrained to 
have zero divergence (see Ref. E] and Section can 
be described by a lattice gauge theory with magnetic 
monopoles-- forbidden. The dipolar correlations can be 
understood by proceeding to a long-wavelength descrip- 
tion with the link variables replaced by a coarse-grained 
vector field.EHU 

Starting from a nearest-neighbor model of spin ice^ 
applying perturbations (for example, magnetic field,'^^!'!^ 
pressurejlSSl impurities ,123 or sim ply the residual further- 
neighbor dipole interactiorpl221) can lead to transitions 
at sufficiently low temperature. The resulting ordered 
phases, which have a finite correlation length and struc- 
ture depending on microscopic details, cannot be de- 
scribed by the coarse-grained theory of the Coulomb 
phase. (The effectively discrete degrees of freedom in 
spin ice should be contrasted with O(A^) spin systems, 
where more conventional critical theories have been de- 
rived from a coarse-gr ained theory.l^^lSil) 

In the present work, these ordering transitions are 
analysed starting from a noncompact U(l) lattice gauge 
theory, in terms of which they are (AndersorP^) Higgs 
transitions involving condensation of an emergent mat- 
ter field (dual to the magnetic monopoles). This per- 
spective leads to a long-wavelength description captur- 
ing both the Coulomb and ordered phases, which plays 



the role of the Ginzburg-Landau theory for conventional 
ordering transitions.^^ This description is strongly con- 
strained by the symmetries, physical and emergent, of 
the gauge theory, leading to nontrivial predictions for 
the allowed ordered states and transitions. 

Of particular theoretical interest are potential contin- 
uous Higgs transitions, which belong to the class of tran- 
sitions beyond the Landau-Ginzburg- Wilson paradigm 
that has recently been the subject of considerable 
attention.'23t2ll Numerical simulations are better suited 
to determining the order of a given transition, but the 
focus here will be on cases where continuous transitions 
appear most plausible. These are argued to include tran- 
sitions into the spin spiral illustrated in Figure [5] and the 
bond-ordered paramagnet of Figure [Sj see Section |Vn| 
for a summary of the transitions. The resulting criti- 
cal behavior may be observable in experiment, but the 
particular forms of the applied perturbations and the si- 
multaneous requirements of low monopole density and 
ergodicity mean that testing these predictions is likely 
more feasible using numerical simulations. 

Previous studies of phase transitions in spin ice have 
considered the low- temperature ordered states favored 
by the dipolar interactions ,'^2] magnetic fields, ^^"^^ and 
pressure.'^S' While such transitions can be included in 
the framework described here, the present work excludes 
states with nonzero net (spin) polarization (see Sec- 
tion [n]); these include those favored by application of 
a uniform magnetic field.'^^'^ Such states require a gen- 
eralization of the method, which will be presented else- 
where. 

The present analysis has connections with studies of 
the cubic dimer modeP^MSJ ^^-^^ quantum models on the 
pyrochlore lattice,'2fi' both of which exhibit Higgs transi- 
tions. An interesting feature of spin ice, distinguishing 
it from these precedents, is symmetry under spin inver- 
sion, to be referred to as "time reversal" . As shown in 
Sections |HI B 2| and |IVB[ this causes extra "Kramers" 
degeneracies, manifested in the hopping model studied 
in Appendix [b] as Dirac cones in the spectrum (see Fig- 
ure [lO| . These imply that the matter modes of relevance 
to the Higgs transitions are more numerous than in cases 
studied previously, leading to a greater diversity of or- 
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dered phases and transitions. 

This work will treat a classical statistical model, ap- 
propriate for spin ice over a range of temperatures, but 
is also of relevance to quantum systems. The long- 
wavelength description of these three-dimensional (3D) 
classical transitions is identical to certain quantum phase 
transitions in 2D systems.!^ (It is in some cases possible 
to make the relationship explicit by p erformin g the map- 
ping in terms of microscopic 

models.EsElSIl) a second 
connection is that 3D quantum models with phase tran- 
sitions at nonzero temperature can be described by clas- 
sical actionsp^ these include transitions in spin ice where 
quantum effects lead to ordered phases.l^Hl The classical 
physics of spin ice is also of relevance to the quantum 
Heisenberg model on the pyrochlore lattice, where ar- 
rangements of singlet dimers can be mapped to spin-ice 
configurations 
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The route from the microscopic model, introduced in 
to a field-theoretical description is chosen to 
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emphasize the nature of the Higgs transitions. The map- 
ping of Section |TT] leads to a lattice gauge theory contain- 
ing a matter field defined on the sites of a diamond lattice, 
and the transitions of interest occur when certain modes 
of this field condense. Because this is a strongly interact- 
ing model, these cannot be determined exactly but can 
be labeled by the irreducible representations of the sym- 
metry group. The symmetries, discussed in Section [TiTl 
are modified in an essential way by the presence of the 
gauge field; the representations themselves are treated in 
Section HVl 

In Section |V) the physical observables of the spin model 
are related to the fields appearing in the gauge theory, 
using the standard procedure of identifying quantities 
with identical behavior under the symmetries. The Higgs 
transition is naturally described in terms of the matter 
fields, and it will be shown that order parameters for cer- 
tain phases can expressed as functions of these. Once a 
physical quantity is associated with each bilinear, one can 
construct an expansion, similar to a Ginzburg-Landau 
theory, in powers of the fields and derivatives, involving 
all terms allowed by symmetry. (In contrast to the stan- 
dard Landau paradigm, the primitive fields that appear 
are not the order parameters for the transition, but in- 
stead emergent matter and gauge fields.) This expansion 



is found in Section VI allowing the ordered states and 
transitions to be determined. A summary of the impor- 
tant results and conclusions is given in Section [VII[ 



B. Spin ice 

"Spin ice"'^ refers to a class of frustrated magnetic 
materials, such as Ho2Ti207 and Dy2Ti207, with lo- 
cal moments S on the sites r of a pyrochlore lattice. 



FIG. 1: The pyrochlore lattice shown as a network of corner- 
sharing tetrahedra (orange and blue), along with the dual di- 
amond lattice (dashed Unes). In spin ice, magnetic moments 
are located at the sites of pyrochlore, which are at the vertices 
of the tetrahedra. A strong easy-axis anisotropy requires the 
spins to align with the axes connecting the centers of neigh- 
boring tetrahedra; dotted (red) lines illustrate these for the 
sites of a single tetrahedron. 



shown in Figure [T] In these materials, the magnetic 
moments are large (^ 10 Bohr magnetons) and expe- 
rience strong local easy-axis crystal fields, so an accu- 
rate microscopic description is provided by a classical 
model of Ising spins with a com bina tion of exchange 
and long-range dipolar interactions.!^ The behavior on 
length scales larger than a few interatomic spacings is 
nonetheless accurately described by a model with purely 
nearest-neighbor ferromagnetic interact ions, as a result of 
a remarkable "projective equivalence" .liflliS Specifically, 
a nearest-neighbor model has an extensively degenerate 
minimal-energy manifold that is identical to the set of 
low-energy states of the full interaction potential, known 
from comparison of experiment with simulations .'^ISHHl 
(The small energy differences between these states in the 
latter case are, along with quantum effects, negligible for 
a range of temperatures.) 

With this is mind, one can write the configuration en- 
ergy as 



H=-J J2 Sr-Sr'+ ViS) 



(1) 



{r.r') 



where J > and the local crystal fields constrain the 
spins to Sr — SrSr, with Sr = ±1- The static unit 
vector Sr lies along the axis joining the centers of the 
two tetrahedra that share the site r, as shown in Fig- 
ure [T] and can be chosen always to point outwards from 
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tetrahedra of a given orientation. The term V{S) is a 
perturbation that may represent coupling to an exter- 
nal field or additional interactions between spins (such 
as those beyond nearest neighbors). Specific examples 
will be given in Section VI these include cases where V 
preserves the full symmetry of the pyrochlore lattice as 
well as the Ising symmetry under Sr —Sr, and those 
where it breaks certain symmetries. Most of the analysis 
will require only that \V\ <C J, where \V\ denotes the 
scale of couplings in V. 

The first term of Eq. ([T]) can be rewritten, up to a 
constant, as 



(2) 



where the outer sum is over all tetrahedra and r G i 
denotes all sites in tetrahedron i. Inspection of a single 
tetrahedron shows that there are 6 configurations for its 
spins that maximize |X) "^i-l, with the vector sum aligned 
with one of the 6 cubic directions. In the pyrochlore 
lattice, neighboring tetrahedra are constrained by their 
shared spin, so the total number of allowed configura- 
tions is less than 6^/^ (for a lattice with N spins), but 
still grows exponentially with N and is in fact closely 
approximated by Pauling's estimate of (3/2)^/^. 

In the present work, it will be assumed throughout 
that the temperature T ^ \V\ ^ J, so the system is con- 
strained to states within this minimal-energy manifold. 
These are said to obey the "ice rule" : a tetrahedron with 
total spin maximal must have two spins pointing inward 
and two pointing outward. In the absence of the pertur- 
bation V, the partition function is given in this limit by 
a sum with equal weight over all ice-rule configurations. 
It is well estab lished, based on theoretical arguments and 
numerics jSluHUl ^j^^t such a sum implies long-range cor- 
relations between the spins of the form 




Sr., 



FIG. 2: Two tetrahedra of the pyrochlore lattice (orange and 
blue), along with the diamond sites (points) situated at the 
centers of each tetrahedron. Each site of pyrochlore, equiva- 
lent to a diamond link, is surrounded by a (nonplanar) hexag- 
onal plaquette of the dual diamond lattice (solid black and 
dashed red lines), so a tetrahedron is enclosed in a "cage" 
consisting of 4 hexagons and 10 diamond sites. The dual- 
diamond links A on which the static gauge field a\ takes the 
value I are shown as dashed (red) lines; the remainder have 
a\ = 0. Every plaquette has an odd number of links with 
a\ ^ and so net flux | (up to an integer), as required by 
Eq. (|8|. (This choice of ax is shown for a larger region of the 
dual diamond lattice in Figure [s]) 



the effective gauge theory, and this perspective allows a 
description that captures both the Coulomb-phase corre- 
lations and the microscopic details that are crucial to an 
ordered state. 
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These are the signature of a "Coulomb phase" , and can 
be found using an effective coarse-grained description in 
terms of a gauge theory.'^l Such correlations have indeed 
been observed in neutron-scattering experiments.l^Hll 

The perturbation V acts to distinguish the ice-rule 
states (i.e., give them different Boltzmann weights), and 
so is effectively projected to this manifold. For \V\ > T 
it can lead to an ordered phase with a finite correlation 
length, and it is the transitions into these phases that 
will be the subject of this work. Although the transition 
may involve the formation of symmetry-breaking order, 
a standard Landau theory written in terms of an order 
parameter is not applicable, since it fails to account for 
the Coulomb correlations on the high-temperature side.!^ 
(Here, "high-temperature" means \V\ < T J; for 
T ^ J the system is a conventional paramagnet.) Such 
transitions can instead be viewed as Higgs transitions in 



II. GAUGE THEORY 

The Coulomb correlation functions of Eq. ^ can be 
derived by describing the spin model in terms of a coarse- 
grained gauge theory.^ In order to preserve those details 
of the microscopies that are crucial to the ordering tran- 
sitions, this section will instead present a mapping to a 
gauge theory on the lattice. As illustrated in Figure [2] 
the sites r of pyrochlore are equivalent to the links of 
a diamond lattice (pyrochlore is the "medial lattice" of 
diamond-). For each directed link £, define Bg taking the 
value ± I if the corresponding spin 5'^ is aligned parallel 
or antiparallel to the link. 

When T ^ J, only states within the ice-rule manifold 
are accessible: the spins are arranged so that two spins 
point into and two out of every tetrahedron of the py- 
rochlore lattice. Such an arrangement obeys the simple 
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condition 



div, B = 0, 



(4) 



where div^ denotes the lattice divergence at the diamond 
site i, given by the sum of Bi on hnks £ pointing out- 
wards from site i. This description in terms of an effective 
magnetic field provides an immediate motivation for the 
gauge-theoretical description, and upon coarse-graining 
implies the correlation function of Eq. ^ . 

Since this work will address the case when the ice rule 
is strictly enforced, Eq. Q will be treated as a constraint, 
and the configuration energy is given simply by 



Hb - V{B) , 



(5) 



where the function V appearing in Eq. ([T]) has been 
rewritten in terms of B. The partition function is given 
by summing over all configurations of B consistent with 
Eq. Q with Boltzmann weight e~^«/^. 

To resolve the constraint of Eq. Q, one can express 
B in terms of a vector potential. The diamond lattice is 
self-dual, so define A on the (directed) links A of a dual 
diamond lattice by 



Bp 



curl^(<>) A , 



(6) 



where curl^(£) is the lattice curl around the plaquette 
(nonplanar hexagon) 7r(£) dual to the link £ (see Fig- 
ure 

Expressing B as in Eq. ^ automatically incorporates 
the constraint on its divergence, but to rewrite the prob- 
lem in terms of A, one must also allow for the restriction 
of Bp to the values ±i on each link £. This will be done 
in a way that may appear unnecessarily complicated, but 
that makes the nature of the Higgs transition directly ap- 
parent. One can allow Ax to take all real values by adding 
constraining terms to the energy. 

Ha = -A^cos27r(AA - ax) + U^(cm:\^ Af 

A IT 

+ V{cnv\A). (7) 

The first term involves the static background a, which is 
defined so that 



curl^ a = 



1 



(8) 



up to an integer, for all dual plaquettes tt. For large A, 
the first term in Eq. ([T]) ensures that Ax differs from ax 
by an integer on every link A, and hence that curler A is 
a half-integer for every plaquette tt. The second term 
then selects the values ±i, so that Ha reproduces the 
original model when the parameters obey A 3> ?7 3> 
|y I , T. The partition function is given by integrating over 
Ax on every link A and is identical (in the appropriate 
limit) to that given by a sum over configurations of B. 

Any choice may be made for the background field 
a, provided that Eq. (|8| is satisfied. This expression 




FIG. 3: The dual diamond lattice, showing a configuration of 
the static gauge field ax (solid black and dashed red lines, as 
in Figure [2| that obeys Eq. Q. A hexagonal plaquette and 
one tetrahedron of the pyrochlore lattice are also shown (in 
yellow and blue respectively); see Figure [l] The Cartesian 
coordinate axes are labeled Ui, and the large points show the 
origin of the coordinate system (black) and the three unit 
vectors e.j (red, green, blue) of face-centered cubic (FCC). 
Diamond consists of this Bravais lattice decorated with a two- 
site unit cell, with sites at and Sq. The large shaded cube 
is the unit cell of the simple cubic lattice corresponding to 
the FCC lattice. Arrows on certain diamond links indicate 



the chains appearing in Eqs. (201 and (231 (single and double 
arrowheads respectively) . 



clearly requires that a be nonuniform, and it in fact 
has important implications for symmetry, as discussed 
in Section 



in 



(There is a close relation with the Hof- 
stadter problem]^ which is explored in more detail in 
Appendix [b|) Figures [2] and [s] show a particular choice 
that obeys Eq. ^ and has the smallest possible unit cell, 
containing two unit cells of the diamond lattice. 

The definition of A in Eq. ^ implies that there is a 



redundancy under 



Ax^ Ax^ gradx X , 



(9) 



where X is defined on the (dual) sites j, and grad is the 
lattice gradient. The constraints in Eq. ([7| imply that 
this gauge symmetry is discrete, given fixed background 
a, despite the continuous nature of A. (The redundancy 
also means that any choice of background obeying Eq. ([s]) 
is equivalent.) 

This gauge symmetry can be made continuous by in- 
troducing new degrees of freedom 9j , defined on the sites 
i of the dual diamond lattice, that are "minimally cou- 



5 



pled" to the gauge field A. The new energy function is 
HA,e = -^J2 cos[27r(ylA -ax)- grad^ 0] 

X 

+ ;7^(curU^)2 + l^(curM), (10) 



states with nonzero overall polarization in any direction; 
perturbations that favor them, such as uniform magnetic 
fields, will be addressed with this approach elsewhere. 



III. SYMMETRIES 



with a continuous gauge symmetry 



A, 



Ax + grad^ X 
6, + 2ttX. 



(11) 
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The latter expression implies that 9 is redundant and its 
introduction does not change the value of the partition 
function. 

Despite this redundancy, it is convenient to include the 
"matter field" 9, giving an XY model coupled to a U(l) 
gauge field.l^ A Higgs transition occurs when T is de- 
creased and the matter field 9 "condenses" , leading to a 
quadratic term in the effective action for the gauge field A 
and the loss of the power-law correlations of the Coulomb 
phase, Eq. ([s]). The background field a, which imple- 
ments certain microscopic details of the original model, 
implies that the condensate configuration is nonuniform. 
The Higgs transition therefore also spontaneously breaks 
the spatial symmetry, and the spatial order is manifest 
in gauge-invariant combinations of the fields 9 and A. 

A description of the Higgs transition is given in terms 
of those modes of the matter field that are most impor- 
tant near the transition. One can eliminate ("integrate 
out") others, renormalizing the parameters appearing in 
the action. This will have the effect of softening the con- 
straints enforced by the parameters A and U in Eq. (10 1, 



without modifying the phase structure or the universal 
behavior near any continuous phase transitions. Simi- 
larly, the matter field will be replaced by a complex vari- 
able ipj ~ e'^^ j whose modulus will not be constrained. 
In these terms, Eq. ( 10 1 is a classical analogue of bosons 



hopping on a lattice; the latter model is studied in the 
Appendix |Bj (It should be noted, however, that the lat- 
tice is three-dimensional, and so this is not the classical- 
quantum mapp ing used previously to study transitions 
in spin ice! ^^ * 



The precise form of the relevant modes cannot be de- 
termined analytically in this strongly coupled system, 
but it is possible to categorize them in terms of repre- 
sentations of the symmetry group. This will involve a 
treatment of the modification of the spatial symmetries 
due to the background field a, with which the next sec- 
tion will be concerned. 

While Eq. ^ can describe any configuration obeying 
the constraint Eq. Q , it is in fact inconvenient in the case 
of a nonzero total spin polarization. It is then preferable 
to replace Eq. ^ hy Bi = bi + curl^(£) A, where 6 is a 
static magnetic neld. The latter quantity is conserved by 
any local dynamics within the ice-rule manifold and acts 
as an additional flux on the matter fields V'j i modifying 
the analysis that follows. The present work will exclude 



As for more conventional Ginzburg-Landau theories, 
symmetries provide an important tool in constructing the 
long-wavelength action. In the present case, the modifi- 
cations due to the background gauge field a are particu- 
larly important. These refiect details of the microscopic 
physics that can be neglected in the coarse-graining pro- 
cedure appropriate to the Coulomb phase, '^^H^^ but are 
crucial for ordering transitions. 



A. Lattice symmetries 



The effective model that will be treated in the follow- 



ing sections is Ha,9, defined in Eq. ( 10 1, whose degrees of 
freedom occupy the sites and links of a diamond lattice. 
The discussion of the spatial symmetries will therefore 
be couched in terms of diamond, but the same symme- 
tries apply to pyrochlore, its medial lattice. (It should 
be recalled from Section |ll] that the diamond lattice on 
which 9 is defined is dual to the original diamond lattice, 
whose links correspond to pyrochlore sites.) 

The Bravais lattice of diamond (see Figure [s]) is face- 
centered cubic (FCC), with primitive vectors ej for i £ 
{1, 2, 3}, related to the Cartesian unit vectors Ui by ei — 
2(1x2 + U3) and cyclic permutations. (In these units, the 
cubic unit cell of FCC has side 4 and the nearest-neighbor 
distance in the pyrochlore lattice is -^2.) Translations 
by the lattice vectors are effected by the operators 
Ti. In the pyrochlore lattice, these translations map a 
tetrahedron onto another of the same orientation. 

The diamond lattice has two sites within the unit cell, 
labeled by sublattice a G {0, 1} and separated by dis- 
placement ^0 = Mi+M2+'it3 ~ |:(ei + e2 + e3). (The two 
sites correspond to tetrahedra of opposite orientation in 
the pyrochlore lattice.) A site on sublattice a has 4 near- 
est neighbors, all on the opposite sublattice, at displace- 
ments (— 1)"^^^, for fj, £ {0, ... ,3}, where 6i = So — Bi. 

The crystallographic space group <S of the diamond 
lattic^^ is denoted O^, while the corresponding point 
group is Ofi- In other words, if the action of each lattice 
symmetry Q on a position vector r is given by r — > Qr = 
MQr -I- Vq, then the rotational parts Mq form the group 
Oh- As for any space group, the pure translations T form 
an invariant subgroup X, and any element of © can be 
constructed as the product TQ of some Tel and one 
of the "coset representatives" Q S 0. The elements of 
Q can be found by associating with each element Mq of 
Oh the appropriate translation Vq, unique up to a lattice 
vector, that maps the lattice onto itself ."^ Because of the 
two-site unit cell, the diamond lattice is nonsymmorphic 
(i.e., Oh ^ O^), having Vq = ^0 for certain elements of 
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Oh- For such transformations, which exchange the two 
diamond sublattices (or two orientations of tetrahedra in 
the pyrochlore lattice), let aq — 1; otherwise ctq = 0. 
The 10 classes of Oh are listed in Table |Tj along with a 
description of the corresponding transformations of the 
diamond lattice. 

Perturbations such as an applied magnetic field or uni- 
axial pressure can reduce the symmetry to a subgroup 
of 25. This will modify the couplings allowed to appear 
in the long-wavelength description and lead to different 
classes of transitions and critical behavior, as discussed 
in Section IVll 



equation Eq. (16) specifies TTQ{j) for all sites j, up to an 



arbitrary uniform phase factor. Note that both ttq and 
uj are phase factors, so 7rQ(j) = 7rg"^(j) and — > j) = 
uj-'^{i ^- j) = uj{j i). 

For the particular case of Eq. ([s]) , it is in fact possible 
to choose the background field so that a\ takes the values 
0, -|-| and — I on every link A, as illustrated in Figure |3] 
In this case uj and ttq are both real, but this assumption 
does not lead to any significant simplification and will not 
be made in the following. (Note that the fact that a real 
gauge exists implies that all gauge-invariant quantities 
are real.) 

One can extend the definition of lo to chains of links. 



B. Magnetic symmetries 

With any choice of the background gauge potential a, 
the fields A and 9 (or ip) obey a modified set of symme- 
tries, which will be referred to as the "magnetic symme- 
try group" (MSG). Emm For each element Q G 0, there is 
a corresponding gauge transformation that must accom- 
pany Q to give a symmetry of Eq. ( 10 ). 

The first term of Ha.b (other terms are trivially sym- 
metric) can be rewritten in terms of ipj 



e'^J as 



A 



,27ri(Ai 



(12) 



where the sum runs over directed links j — >■ j of the dual 
diamond lattice, with each link counted twice. Under a 
symmetry Q that maps the site j to Qj, the degrees of 
freedom transform according to 



SO Eq. ( 12 ) becomes 



E 



(13) 



(14) 



after a change of the summation variables. 

The transformation Q is therefore not by itself a sym- 
metry of H^^\ but must be accompanied by a gauge 
transformation. For every Q G 0, define a corresponding 
operator Q under which 



A, 



— An- 



(15) 



where the phase factor ttq compensates the change in 
background field induced by Q. Comparison of Eqs. ( 12 ) 
and ( 14 1 shows that one must choose ttq so that 



^Q(j>Q(i) = ^ jXQ* 



(16) 



2'jriax 



where i and j are neighboring sites and w(A) = e 
Given a fixed choice of background field, the defining 



uj{i j 



£) = uj{i ^ j)uj{j 



allowing Eq. (161 to be generalized to 



, (18) 



where the sites i and j need not be neighbors. The sites 
visited by the two chains in Eq. ( 18 ) must be in corre- 



spondence, but the path is otherwise arbitrary. 

It should be noted that uj{i —>■•••—>■ j) depends on the 
route taken by the chain, because the background field 
a has nonzero lattice curl. In fact, because curia = ±| 
for all hexagonal plaquettes of the dual diamond lattice, 
any closed path has uj{i — > • • • — > i) = (—1)^'" where 
iVh is the number of plaquettes it encloses. This gauge- 
independent result allows the commutation relations of 
the operators Q to be found. 

Consider for example the elementary translation oper- 
ators Ti, which obey T2T1 = T1T2. To find the commu- 
tation relations of the corresponding MSG operators %, 



consider tt. 
relation 



T — — T 
'2 '1 ' 2 ' 1 



(j). The definition of ttq gives the 



(19) 



which, together with Eq. (18), gives 



(j) = io{j ^ Tij ^ T1T2J ^ T2J 



■3). 
(20) 

using the convention that, for second neighbors i and j, 
i ^ j denotes the chain through their common neighbor. 
This expression involves a closed chain and so is gauge 
invariant; inspection of the example in Figure |3] shows 
that it encloses a single diamond plaquette. This is in 
fact true for j on either diamond sublattice, so one has 
t:^ — — ij- J- (j) = —1 for any site j, independent of the 
gauge.' 

This is an instance of the general result that the MSG 
operator Q corresponding to a space-group operator Q is 
fixed by Eq. ( 16 ) only up to the arbitrary global phase 



of TTQ. (As a special case, for any MSG operator £ corre- 
sponding to the identity transformation E, TTs{j) is equal 
to a uniform phase, independent of j.) If — Q denotes the 
MSG operator corresponding to the same space-group el- 
ement as Q but with 7r_Q = — ttq, then one can write 



T2T1 



-T1T2 



(21) 
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TABLE I: Classes of the point group Oh,*"^ and corresponding symmetries Q of the diamond and pyrochlore lattices. For each 
class, det Q = ±1 for proper and improper transformations respectively, and ctq = 1 for those that exchange diamond sublattices 
(or pyrochlore tetrahedra of opposite orientations). In each case an example element is listed using "XYZ" notation:^ the 
order of the letters indicates the permutation given by the matrix Mq and bars signify minus signs. 



Class 


Order 


detQ 


CQ 


Example 


Description 


E 


1 


1 





XYZ 


Identity (E) 


SCs 


8 


1 





YZX 


Proper rotation (R) by ±4^ about diamond link Si 


3C2 


3 


1 





XYZ 


Proper rotation by tt about cubic axis Ui 


6C4 


6 


1 


1 


XZY 


Screw rotation by ±^ about cubic axis Ui 


6^ 


6 


1 


1 


XZY 


Screw rotation by n about FCC vector 


I 


1 


-1 


1 


XYZ 


Inversion 


8Se 


8 


-1 


1 


YZX 


Improper rotation by ±^ about diamond link Si 




3 


-1 


1 


XYZ 


Glide reflection in plane perpendicular to Ui 


65*4 


6 


-1 





XZY 


Improper rotation by ±^ about cubic axis Ui 


6crd 


6 


-1 





XZY 


Reflection in plane normal to FCC vector 



One similarly finds T1T3 — —T3T1 and 7372 = —T2Tz- 
The commutation relations of these translations with 
other symmetry operations are calculated in AppendixjX] 



1. Cubic translations 

Representations of the MSG are most easily found us- 
ing an Abelian invariant subgroup of translation opera- 
tors. The primitive FCC translations Ti do not commute 
with each other, so consider instead 



(22) 



This combination ensures that [/Ci,7}] = for all i and 
j, and so [Ki^fCj] = 0. The corresponding space-group 



operators K^, which will be referred to as "cubic trans- 
lations" , translate by Aui , the unit vectors of the simple 
cubic lattice upon which FCC is based. (The cubic unit 
cell, illustrated in Figure |3j contains four FCC sites and 
hence eight sites of the diamond lattice.) 

It is shown in Appendix |A] that one can choose the 
global phases so that 



(23) 



where the chain j — s- K^j follows a right-handed helix; 
see Figure [3] and Eq. (A13|. This allows the commuta- 
tion relations of K-i with other any MSG operator Q to 
be found. Suppose some space-group operator Q obeys 
QKiQ^i = K'; K' is then a cubic translation by 4MQMi. 
The commutation relation can be inferred from 



= 7r^'(j)a;*(j ^ • •• - 



(24) 
(25) 



By choosing a right-handed helix for the first chain, which 
spans a cubic translation from j to K'j, one can make 
use of Eq. (23). For a proper transformation Q, the sec- 



ond chain is also then a right-handed helix and one has 
^yc'-iQKiiQ-i (i) = 1- By contrast, an improper transfor- 
mation reverses chiralities, and inspection of the diamond 
lattice shows that the right- and left-handed spirals differ 
by their path around a single hexagonal plaquctte. This 
leads finally to the important conclusion that 



Q/C, =/C'Q det ( 
where detO = det Mo = ±1. 



(26) 



2. Time-reversal symmetry 



Besides the geometrical transformations of the py- 
rochlore (or diamond) lattice, another important symme- 
try of spin ice is that under reversal of all spins, which 
will be referred to as "time reversal" and denoted 9. As 
remarked above, this symmetry will be of profound im- 
portance for the eventual conclusions. In terms of the 

variables Bg, Q obviously amounts to Bi i Bi (for all 

diamond links i), which preserves the spin-ice constraint 
of Eq. ^ and is a symmetry of the configuration energy 
Hg provided V contains only even powers of B. (This 
will not be the case if, for example, an external magnetic 
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field is applied to the spins.) 

Eq. ^ gives, correspondingly, Ax ^ — ^a, under 
which the last two terms of H^fi, in Eq. (10), are again 
trivially symmetric (for ap propriate V). Considering 



once more H'-^^ in Eq. (12), it is necessary to transform 
ipj according to 



(27) 



so Q is an antiunitary operator. Time reversal is its own 
inverse, so 6^ = 1, implying that |7re(j)| = 1. The phase 
should be chosen such that 

TTe (z)^e 0>* (* ^ i) = ^* (j ^ > (28) 

and the terms in H^-^^ for i ~^ j and j ^ i are exchanged 
by &. Note that if the background gauge potential is 
chosen real, then uj{i — > j) = oj{j — > i) and one can set 
TTe = 1. 

Applying 9 and the MSG operator Q successively to 
ipj shows that the two commute only if 

7rQ(Q~^J>e(Q~V) = vre(j)7r^(Q"ij) , (29) 
or equivalently, using Eq. ([28|, if 



1 . 



(30) 



for all sites j. In fact, it is straightforward to show, using 
Eq. (18) and the fact that uj{i ■ ■ ■ i) — ±1, that 
if this is true for any site, it is true for all sites. By 
choosing the global phase of ttq so that 7rQ(0)a;(0 -> QO) 
is real, one can therefore ensure that Q commutes with 
any MSG operator Q. This condition is indeed satisfied 
by the conventions used here for the translation operators 
Ti and ICi, according to Eqs. (A8) and (A15). 



IV. MSG REPRESENTATIONS 

The Higgs transitions involve condensation of the mat- 
ter field Ipj, and the nature of the ordered phase is de- 
termined by the configuration of the condensate. While 
the details will depend on microscopies, the symmetry 
structure can be elucidated by considering the irreducible 
representations, or "irreps" , of the MSG, which label the 
modes of the matter field. The irreps will be found by 
following the general procedure for space groups^''- but 
there will be certain differen ces because of the commu- 
tation relations such as Eq. (21).!^ No attempt will be 



made to find all representations systematically; interest 
will be restricted to those with smallest dimension for 
certain high-symmetry wavevectors. 

As remarked in Section [III B 1[ it is convenient to con- 
struct representations starting from an invariant sub- 
group of mutually commuting translation operators, such 
as the cubic translations ICi. The vectors belonging to a 
given representation will be labeled by their eigenvalue 
of the operators ICi, according to 



note that translates by 414^. Distinct representations 
are associated with "wavevectors" k within the region 
where — ^ < k ■ Ui < j for all i. This cubic "reduced 
Brillouin zone" Q3r has a volume of 1/4 that of the Bril- 
louin zone *Bl of the FCC Bravais lattice (a truncated 
octahedron). The notation [k] will be used to signify 
wavevector n reduced to *8r by addition of multiples of 
the reciprocal lattice vectors ^Ui. 



A. FCC translations 

The FCC translations Ti commute with ICi, so one can 
also label the vectors by their eigenvalue under Ti, say. 



From Eq. (22), one has Tf — — /C1/C2, so for given k there 



are only two possible eigenvalues of Ts , 

T3\K,£)=ii-lYe~'^--^\K,£), (32) 

where £ S {0,1}. 

Since the remaining FCC translation operators 7i and 
T2 commute with ICi but not T3, they mix vectors with 
equal k but distinct £. For instance, T3T1 = — TiTs, 
so the vector 711^,^?) has Ts-eigenvalue — i(— l)''c^"' '^^, 
allowing one to define 



which implies 



=_ie---iri|K,0), 



\k,0) = -ie-'"-"iri|/«,l) 



(33) 



(34) 



The action of T2 on the eigenvectors can then be found 
by using 7^ = —IC^T^^T'^ ■ These results can be sum- 
marized as 



7^1/^,^) =i^(7^,,e-"^-^'|K,0, 



(35) 



where a^'^'^ are the Pauli matrices. 



B. Time reversal 

The time-reversal operator O, introduced in Sec- 
tion [TirBj2] is antiunitary and commutes with the trans- 
lation operators % and ICi. Applying O to Eq. (31) one 



therefore finds that 8|k) is also an eigenvector of ICi, but 
with wavevector [—/«]. (Reduction to *Br is only strictly 
necessary for k at the edge of *Br, but these points will 
be of particular importance.) 

One can therefore define the matrix 6{k), 



e|/«,£) = ^0r£(K)|[-/^],n 



(36) 



The commutation of and % then gives 



ICrln) 



g-i«.(4«,)|;^\ 



(31) 



cr'6»(K)e~ 



-6»(K^)(cr')*e" 



(37) 
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for all i G {1,2,3}. This is sufficient to fix 9{k) up to 
phase, but 6{k) in fact depends on the relation between 
[— /«] and /t, and can be written as 





if 


— K 


7 2 

cr-'cr 


if 


— K 




if 


— K 


,t2 


if 


— K 



= —K — 



K (body of *Br) 
K— ^Uj (face) 
fcj (edge) 
— |(5o (corner), 
(38) 

1, which relates 



where |^9('«)| = 1- 

For consistency, one requires 
the values of 'd{K) and i}(—k) at a generic wavevector k. 
For the points k — and n = — f^o (the center and 
corner of *Br), however, Eq. (38) gives a contradiction. 
This implies that there must be an enlarged (Kramers) 
degeneracy at these points, and another label besides i 
is required. The hopping model studied in Appendix [B] 
has Dirac cones in the spectrum around these points, 
confirming this result. 



C. Other transformations 

The action of the remaining MSG operators Q on a 
vector \k,£) can be determined using the commutation 
relations of Q with % and ICi . Treating first the k label, 
acting with Q on the defining equation Eq. ( 31 1 and using 
Eq. (|26|) gives 



JC'Q\k,£} det( 



-ire-(4ui) 



(39) 



The operator /C' translates by AlsAqUi, so Q\k,£) 
has wavevector /«q(k), where e""*^^'^-''^^''^^"'-' detQ = 
g-ire-(4iii)^ For detQ = 1, the wavevector therefore trans- 
forms as usual under the space-group operation Q, with 
kq{k,) — [Mq/«]. For improper transformations, the fac- 
tor of det Q = — 1 appearing in the commutation relation 
leads to /«q(«) = [Mq/« — jSq]- 

The set of distinct wavevectors given by kq{k) for all 
Q is referred to as the "star" of k. Labeling the 
wavevectors within a star by for k G {1, 2, . . . , |*t|}, 
one can write 



(40) 



defining the matrix A, with 2|7k-K| rows and columns. Ap- 
plying a second MSG operator Q' to Eq. (40 1 gives 



A(Q'Q) = A(Q')A(Q) 



(41) 



so these matrices provide a unitary representation of the 
MSG. (This is also therefore a projective representation 
of the space group, with exactly the same factor system 
as the MSG.) According to Eqs. (31) and (35), one has 
A«'^',.£(/CO =e-''^-(4«')(5««.(5«. and 



it is similarly useful to define A(0). 

The matrix A(Q) is in principle determined by the 
commutation relations of Q with Ti, but this calcula- 
tion is complicated by the need to make a consistent as- 
signment of the relative phase of \k,£) between different 
wavevectors k^. It is simpler to use the transformation 
properties of the eigenstates of the hopping model treated 
in Appendix [B] which has the same symmetry group and 
so the same set of representations. The main interest, 
however, is not the transformation properties of the vec- 
tors \k,,(), but rather those of the gauge- invariant combi- 
nations of these, to be introduced in Section [V A| These 
can in certain cases (the diagonal bilinears) be calculated 
using only the commutation relations of the matrices A, 
which are identical to those of the MSG operators Q. 



D. Relevant matter-field modes 

The stars */« play an important role in the long- 
wavelength description: if a mode of the matter field 
with wavevector k is relevant at a transition, then sym- 
metry requires that the corresponding mode at /^©(h;) is 
also. The set of matter fields appearing in the action 
will therefore generically consist of a single star ★/«. If 
the symmetry group is reduced by the perturbations V 
in Eq. ([T]) , then the number of wavevectors in a star can 



A 



(42) 



be correspondingly reduced (see Section VIA). 

For a generic point, -kn has 96 members (the point 
group has 48 elements; time reversal contributes a factor 
of 2), but certain high-symmetry points have consider- 
ably smaller stars. The points k — and —jSq are linked 
by improper transformations and form a star consisting 
of only two wavevectors. As noted in Section pVB[ how- 
ever, time-reversal symmetry increases the dimensional- 
ity of the corresponding representations and the pres- 
ence of Dirac cones at these points in the hopping model 
studied in Appendix [B| makes Higgs transitions involving 
these modes appear unlikely. 

The only other wavevectors that are invariant under 
time reversal are — and — fe^, i G {1,2,3}, which 
form separate stars with respect to the space group but 
are linked by improper MSG transformations. (The two 
sets are equivalent and can be exchanged by choosing 
the opposite sign in Eq. ([22]).) They therefore form a 
star with 6 members, to be denoted -^Kq; all other stars 
consist of more than 6 wavevectors. 

As described in detail in Section |V) the symmetry of 
an ordered Higgs phase depends on that of the wavevec- 
tor(s) K at which condensation occurs. Of most interest 
in spin ice are those ordered states with definite sign un- 
der time reversal (odd for spin-ordered, even for param- 
agnetic); such states require that condensation occurs at 
time-reversal-symmetric wavevectors. The concrete re- 
sults presented here are therefore based on the assump- 
tion that */«6 is the unique star that is relevant at the 
transition. (These modes indeed have minimal energy 
for a range of parameters in the hopping model of Ap- 
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pendix|B]) Exactly the same procedure can be applied 
when other stars are relevant. 

It is convenient to introduce an alternative labeling for 
the wavevectors in ★Kg: let = ^li'^j ^'^'i ^j- = 
— for j G {1,2,3}. The reason for this notation is 
that, as described in SectionjV] the matter fields at these 
six wavevectors can be used to construct order parame- 
ters for six spin spirals (see Figurejs]), parallel to the three 
cubic axes and with right- and left-handed chiralities. 



so Dfj^ip^i is gauge invariant (but transforms like ip under 
global phase rotations). The operator D is also a vector 
under the space group. 

Under time reversal, the matter-field configuration {ip) 
is replaced by so the continuu m field transforms 

as <^ A{Q)ip*. As noted in Section ; 
under time reversal; the same is true for a. 



Ax is odd 



so 



(48) 



E. Continuum fields 



ORDER PARAMETERS 



The physics of the Coulomb phase is described by a 
long- wavelength theory found by coarse-graining^ the 
gauge field A\. A theory of the Higgs transitions must 
also involve the matter fields V'j, whose mode structure 
carries information about the microscopic details that is 
important in the ordered phase. 

The set of matter-field modes that are relevant to the 
Higgs transition form a representation of the MSG, so the 
notation \ki^,£) will henceforth be used to denote such a 
mode. (In the generic case, no further indices are re- 
quired to label the complete set of relevant modes.) The 
important configurations of the matter field are those 
that are "close to" these modes, and can be written us- 
ing the same notation as 



The continuum fields ¥'(t) and a.{r) provide the de- 
sired long-wavelength description, but the strongly cou- 
pled nature of the system prevents a direct derivation 
of the relation between these and physical observables. 
Each observable can instead be associated with a func- 
tion of the continuum fields with identical symmetries, 
and the two are, to leading order, proportional. 

It might appear most natural to express the order pa- 
rameters in terms of the gauge field A, related to the 
spins S{r) by Eq. ([6|, or its continuum equivalent a. 
The Higgs transition is, however, more readily described 
in terms of the matter fields ip, and it will be shown that 
the order parameters can also be expressed as combina- 
tions of these. Because of the gauge symmetry under 
Eq. (46), the lowest-order combinations that can corre- 



(43) spond to physical quantities are bilinears of ip. 



where Vr is a projector that restricts to a large region 
surrounding the point r. This expression defines the con- 
tinuum field ip{r), a 2|*/«|-component vector, which is 
assumed to be slowly varying on the lattice scale. Under 
a spatial symmetry Q, the configuration is transformed 
to Q\i}}), giving the corresponding transformation 



V{r) ^ A(Q)(p(Q-ir) 



(44) 



for the continuum field. 

The continuum version of the gauge field A\ will be 
denoted a.{r) and transforms as a vector. 



Q.{r) 



(45) 



under the space-group operation Q. The microscopic 
gauge symmetry Eq. ( 11 ) then transforms the continuum 
fields a and ip as 



a a + dX(r) 

ip ^ ipc'^^'-') 



(46) 



(where a factor of 2tt has been eliminated by the defi- 
nition of a). It is therefore convenient to introduce the 
"covariant derivative" 



(47) 



A. Gauge-invariant bilinears 

Given the set of 2|*/«| matter fields (p,^i{r), one can 
construct all bilinears as 



(r) = ipl{r)(T''ip^,{r) 



(49) 



where /i € {0, . . . , 3} and cr° is the 2x2 identity. (In 
this expression, the matrix notation applies only to the £ 
indices.) These quantities are not independent, satisfying 



E$* , $ , = 2$^ <I>n / 



(50) 



From Eq. (44), the bilinears transform according to 

*(r) 4^ A(Q)*(Q- V) , (51) 
where the matrix A is given by 

A,.^.„.«;.^(Q) = lTrAl^,{Q)cT^A^^^,iQ)a'^. 

(52) 

Note that the overall phase of the MSG operator Q has 
no effect on A(Q), and that these matrices form a uni- 
tary (vector, rather than projective) representation of the 
space group. 
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One can therefore associate with each bihnear a 
wavevector referred to the fuU FCC Brillouin zone *Bl- 



Using Eq. ( 42 ) , one finds 



A 



1 ^K.-l 



(53) 

so the bilinear $pkik2 has wavevector ~ '^ki — f 
(where Mq = 0). For a given star ★/«, the *8L-wavevectors 
appearing in the representation provided by A(Q) are 
given by all such combinations of the *8R-wavevectors in 
★/«. These will comprise more than one star under Q (i.e., 
more than one closed set under the action of the matri- 
ces Mq) and so the representation A(Q) is necessarily 
reducible, even when A(Q) is irreducible. 

Of particular interest are the "diagonal" bilinears of 
the form $^kk, where the two wavevectors are identical. 



-1 flKK, 



*o«k: allowing 



These are real and obey 

each bilinear to be treated as a 3-component vector 
with magnitude given by = ^Okk- The wavevec- 

tor associated with $^kk is — f m^, irrespective of k, and 
these are clearly the only bilinears whose wavevector van- 
ishes on reduction to QSr. The diagonal bilinears there- 
fore correspond to order parameters whose unit cell is no 
larger than the cubic unit cell of the pyrochlore lattice 
(containing 16 sites). 

The wavevectors — f are closed under the action of 
the space-group symmetries, so the transformation ma- 
trices A^kk,i/k'k' for this subset form a representation of 
the space group (a "subrepresentation" of that provided 
by A). As noted in Section IV C[ these matrix elements 
can be calculated directly based on the multiplication 



table of A(Q). From Eq. (52) 



X Tr AL,(Q)A_(r^)A„„,(Q)A,,,,(r.) , (54) 

where the Pauli matrices have been expressed in terms 
of A(7i) using Eq. (42). The matrices A form a unitary 
representation of the MSG, so this can be rewritten as 



X TrA,,^,iQ-%Q%) . (55) 

The operator Q-^T^Q% gives a pure translation by 
+ MQe^, and can be simplified using Eq. (All ). The 
matrix element A^ kk.i ^k'k' (Q) is finally found by once 
more applying Eq. (42). 



The transformation properties of $(r) under time re- 
versal can be found straightforwardly using the results 
of Section [IV B[ For the diagonal bilinears involving -kKg, 
for example, one has 



B. Observables 

To identify the bilinears discussed above with appropri- 
ate physical quantities, it is necessary to find an observ- 
able that transforms under the symmetries of the prob- 
lem in the same way as a given bilinear. The first step 
is therefore to identify the potential order parameters in 
the physical model and determine their transformation 
properties. 

The observables of primary interest are of course the 
spin degrees of freedom Sr- These are odd under time re- 
versal (by definition of the latter) and are pseudovectors: 
under a space-group symmetry Q, 



(57) 



One can alternatively make the Ising nature of the spins 
explicit using the scalar variable Sr defined in Section [lB[ 
which transforms as Sr <— (det Q)(— l)'^*SQ-i.r- 

Those bilinears that are even under time reversal can- 
not correspond to spin degrees of freedom, so one is led 
to consider also correlations of the spins. (Such observ- 
ables can be order parameters for paramagnetic ordered 

In particular, the 
±iis 



phases, as discussed in Section VIE 



product of two nearest-neighbor spins ■ Sr 
time-reversal even and transforms as a scalar under the 
space group. Within the ice-rule manifold, every tetra- 
hedron has two edges where this is negative, on opposite 
sides of the tetrahedron. (These are the frustrated fer- 
romagnetic bonds, linking sites with the same value of 
the Ising variable s.) One can therefore combine these 
nearest-neighbor correlations into a nematic variable Ni 
for tetrahedron i, which lies along the cubic direction 
joining the edges with negative correlation and trans- 
forms as Ni ^ zbMQiVQi. 

As illustrated in Figure |4j each site of the pyrochlore 
lattice has 12 second-nearest neighbors, separated by dis- 
placements of ^(ei +€2) and equivalent vectors. General 
combinations of correlations between second-neighbor 
spins will not be treated, but it is worth observing the 
pairs can be divided into two "pinwheels" of opposite 
chirality. While the symmetric combination of all cor- 
relations is clearly a scalar quantity, the difference be- 
tween correlations on right- and left-handed pinwheels is 
a pseudoscalar (changing sign under improper transfor- 
mations) . 

Further-neighbor pairs of spins and higher-order cor- 
relations can also be considered, but these are subject 
to rapidly diminishing returns. One is also limited by 
the restriction to the ice-rule manifold; if this were re- 
laxed, the monopole density would, for example, provide 
a scalar order parameter. 



(-i)'"<i>M,.±..± 



(56) 



where the notation introduced in Section IIVDI is used. 
There is therefore a single diagonal bilinear that is odd 
under time reversal at each wavevector Kj±. 



C. Irreducible representations 

The observables can be combined into a vector F, so 
that any linear combination can be expressed as v ■ F, 
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FIG. 4: The 12 second-nearest neighbors of a site of the py- 
rochlore lattice. The second-neighbor pairs can be separated 
into two "pinwheels" of opposite chirality, joined with solid 
(green) and dashed (red) lines. (The chirality is defined by 
the scalar triple product of the three vectors illustrated with 
arrows: the two solid arrows link a pair via their common 
neighbor and the dashed arrow links the centers of their tetra- 
hedra.) The difference between the spin correlations on the 
two pinwheels gives a pseudoscalar order parameter, labeled 
"-" in Table HU 



where the vector v transforms according to matrices 
D(Q) determined by the considerations of the previous 
section. One therefore has two representations of the 
space group, given by the transformation matrices for 
the bilinears A(Q) and those for the observables. These 
two sets of quantities will be identified by exp ressing both 
in terms of the irreps of the space group.BS 

An irrep T''^ is labeled by a *BL-wavevector k and 
an index p, and consists of a set of dkp x dkp matrices 
r'°^(Q) that obey the multiplication table of the space 
group. Since the translation operators do not com- 
mute with the other operators Q, a given irrep involves 
all wavevectors in the star irk. (The symbol k will be 
used for wavevectors referred to *Bl and -kk for the set of 
wavevectors resulting from the action of Mq on fe.) The 
matrices for translation operators T''^{T) are diagonal, 
with the element in row a given by e~''°" , where ka G -kk. 

The projection matrices given bjEH 



D 



W\ 



K 



D( 



(58) 



a € {1,2,..., dkp} form a set that transform under 
according to the matrix r'^^(Q): 



D( 



b 



)Vb ■ 



(59) 



The observable F^p — Va ■ F is then said to transform 
as row a of the representation T'^^. (While this state- 
ment is formally true for any choice of starting vector v, 
one can certainly choose it such that Va = for all a, 
in which case the result is of course trivial. If there are 
different starting vectors v and v' that give linearly inde- 
pendent results and v'^ , these are said to provide mul- 
tiple "copies" of the representation.) The corresponding 
combination of the bilinears, denoted ^a^, transforms in 
exactly the same way as F^^^, allowing the pair to be 
identified up to a constant of proportionality. 

Time-reversal symmetry can be included by doubling 
the number of representations, with each labeled either 
even or odd under &. (The full group is a direct prod- 
uct of the time-reversal group and the space group; the 
former is Abelian.) Each projection operator will then 
include a factor projecting into the subspace with appro- 
priate sign under time reversal. 



Including this factor and rewriting the sum in Eq. ( 58 ) 



in terms of the subgroup of translations T and the coset 
representatives (defined in Section III A I gives 



D 



kpzk 



1 



, --p(E)±D(e)] 



*D( 



(60) 



where ± denotes an irrep that is even or odd under time 
reversal. The second factor in this product involves an 
infinite sum over translations by all FCC-lattice vec- 
tors r, and has the effect of projecting into the subspace 
of observables with wavevector fc^- 

There are 24 diagonal bilinears that can be constructed 
using the wavevectors comprising ^Kq, defined in Sec- 
tion |IVD[ and the linear combination for each irrep is 
found using the analogue of Eq. (60). Results are pre- 



sented in Table |II| for the diagonal bilinears, which, as 



can be used to find a set of observables that transform 
according the irrep F*^^. (The equivalent set of matri- 
ces for the bilinears will be denoted ^ab\ In particu- Tij/ (in the naming convention of Bouckaertpl 
lar, for any choice of v, the vectors Va = D^fD^^ff for tation Tig of the point group Oh), and the corresponding 



noted in Section [V A[ have wavevectors k — —^u^j^ for 
fi £ {0, . . . , 3}. I'hese can be divided into two stars 
under the space group, to be denoted ★fcr — {0} and 
*kx = { — ^Ui}, corresponding to the symmetry points 
r and X respectively. Observables with these wavevec- 
tors have period no larger than the cubic unit cell of the 
pyrochlore lattice and will provide the order parameters 
for the phases of interest in Section |VI[ 

Linear combinations of the spins transform according 
to representations r'^^" , and there is in fact only one rel- 
evant representation for each of the stars *fcr and ★fex- 
For the former, this is the pseudovector representation 



represen- 
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TABLE II: Order parameters for transitions in spin ice and the corresponding diagonal bilinears constructed from the modes of 
the matter field at wavevectors -kKa = {t^x} where j € {1,2,3} and x = il- The irreps of are labeled using the standard 
Bouckaert convention,*^'- with the superscript specifying the sign under time reversal. The nematic observables, illustrated in 
Figure [g] are labeled djr and djT, where r denotes the structure on the two tetrahedron orientations (correlations on tetrahedra 
of a single orientation for r — 0,1; symmetric or antisymmetric combinations for r = ±). The vectors v are defined by 
v{di) — :^(1, — |, ~|)) v{di) — |(0, 1, —1). All sums are over j £ {1, 2, 3} and x = ij the diagonal bilinear ^i.jxJx i^ denoted 
$i,jx brevity. 



irrep F 

xj 



observables 
+ 

di+, di+ 
di_, di_ 
rfjo(-fnj), dji(-fuj) 

M 

^^7X 



bilinears 

Ejx "j('^i)x*o,jx. Ejx *'j('^i)x*Oj, 



"In the subscript, [j it 1] indicates arithmetic modulo 3. 
''The uniform magnetization does not correspond to any combi- 
nation of these bilinears. 



description 
identity 
pseudoscalar (Figure |4| 



see Figure [5] 

) ^(^j,[J+i]+-^j.lj-i]-T seeFigure|8| 



magnetization 
spin spiral (Figure |5 



observable is the uniform magnetization M . As noted in 
Section [nj the present analysis is restricted to situations 
where there is no net polarization, and so these order 
parameters vanish. 

At the star *fcx, the irrep X3 is associated with spin 
spirals, as illustrated in Figure [s] (Such states were 
found in the Monte Carlo studies of Melko et al.^ 
who described them as consisting of antiferromagneti- 
cally stacked planes of tetrahedra; see Section VII ) This 



representation is 6-dimensional, and the 6 observables, 
denoted Sj±^ are order parameters for right- and left- 
handed spirals along the 3 cubic directions. The corre- 
sponding bilinears are given by simply ^jj±.j± ■ As noted 
in Section [IV C[ improper transformations exchange mo- 
menta K,j± and hence the spiral chiralities. 

The observables for irreps r'^^"'' are somewhat more di- 
verse. Consider first the 1-dimensional representations at 
*kr, scalar Fi and pseudoscalar Fi'. The former is fully 
symmetric and is given (for any star -kn) by the symmet- 
ric combination of bilinears ^qkk- As noted in Sec- 
tion VB an observable transforming as a pseudoscalar 
(invariant under proper transformations, but odd under 
spatial inversion) can be constructed using correlations 
of second-neighbor spins. The bilinear for the same irrep 
can be written X^jxX^Ojxjxi involves the difference 
between the bilinears for the wavevectors Kj^ associated 
with opposite chirality x. 

The remaining observables can be expressed in terms 
of the nematic variable ATj, defined in Section [VB[ mea- 
suring the nearest-neighbor correlations on tetrahedron 
i. For the star *fcr, the relevant irreps are F12 and 
Fi2' (representations Eg and Eu of Oh), which are both 
2-dimensional. The observables corresponding to both 
have spatial structures, illustrated in Figure [6| related to 




FIG. 5: A spin spiral maximizing ei+ and an example of a 
perturbation that favors it. There are 12 spiral-ordered states, 
with either sign for one of the 6 order parameters Sjx describ- 
ing spirals of chirality x = i along the cubic axis Uj. The 
observable labeled di_ in Table [nj is shown by solid (black) 
and dashed (red) lines, indicating respectively positive or neg- 
ative signs for the correlations on these links. This spin spi- 
ral (or its time-reversed partner) also maximizes di_: spins 
joined by solid and dashed bonds have, respectively, positive 
and negative inner product. An additional nearest-neighbor 
Heisenberg interaction with this symmetry (ferromagnetic on 
solid bonds, antiferromagnetic on dashed) therefore has en- 
ergy minimized by the spin configuration shown. 
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FIG. 6: Nematic order parameters labeled dj and dj (where 
the cubic unit vector Uj is vertical) for a single tetrahedron. 
These observables involve the illustrated linear combinations 
of nearest-neighbor spin correlations, where the line thickness 
indicates magnitude and dashed (red) lines have negative sign. 
The symmetries of these order parameters are those of the Eg 
cubic harmonics, d^^ and d^2_y2 (where z is vertical), which 
are plotted at the centers of the tetrahedra. 

the Eg cubic harmonics c?22 and dx2_y2. They will be 
denoted di^ and di^ for the irrep ri2, and di_ and di_ 
(see Figurejs]) for T12'] the subscript ± indicates that the 
sign for the correlations is the same or opposite on the 
two tetrahedron orientations. (As for the cubic harmon- 
ics, a different choice of "quantization axis" can be made 
by taking different linear combinations.) The observables 
for irreps and have wavevector k — —-^Ui and 
correlations restricted to tetrahedra of a single orienta- 
tion. They have a similar c?-wave structure, but with the 
symmetry axis along Ui. 

As is evident from the case of the uniform magnetiza- 
tion, there is in general no guarantee that the observables 
and bilinears be in one-to-one correspondence. The in- 
tention is to include enough observables that there is at 
least one for every bilinear, and it is obviously permissi- 
ble to ignore those observables that do not correspond to 
any bilinear. There might, however, be multiple sets of 
observables transforming in the same way, in which case 
any linear combination can be chosen, since all behave 
similarly. Alternatively, multiple combinations of bilin- 
ears may have the same transformation properties, and 
the observable will be given by some unknown combina- 
tion of the given bilinears. 

VI. PHASE TRANSITIONS 

The nearest-neighbor model remains in a Coulomb 
phase at arbitrarily low temperatures, but perturbations 
can drive it into a Higgs phase. A standard Landau the- 
ory expressed directly in terms of the order parameters 
does not describe such a transition, but one can instead 
write down an effective action in terms of the coarse- 
grained fields ip and a. This contains, in principle, all 
terms compatible with symmetry, and can be expressed 



as a power series in successive powers of the fields and 
derivatives. The terms of most interest will in fact be 
those containing only the fields cp, and these are conve- 
niently expressed in terms of the order parameters 

This relies on the assumption that there is some non- 
trivial long-wavelength description, in terms of smooth 
variations on the minimal-energy configurations. The fo- 
cus in the following will therefore be on cases where this 
is true, viz. for continuous transitions and also for weak 
first-order transitions, provided the correlation length be- 
comes sufficiently large. Like a standard Landau theory, 
such an expansion can nonetheless be useful even when 
the transition is strongly first-order, where it can often 
at least exclude the possibility of a continuous transition. 

Given the long-wavelength action, one can determine 
the possible transitions based on the behavior of the or- 
der parameters as the coupling constants are varied. Al- 
though certain ordered states and the perturbations that 
select them will be described, this should not be consid- 
ered the main focus of the current work. Numerical simu- 
lations are better suited to address many such questions, 
which can be difficult to answer definitively using ana- 
lytic arguments. Where the applied perturbations pre- 
serve the full space group, the approach will be to give 
an outline of the general types of transitions that are pos- 
sible. Then specific examples will be given of symmetry- 
breaking perturbations for which the resulting ordered 
states and phase transitions can be determined. 



A. Long- wavelength action 

The aim is therefore to find a long-wavelength action 
(free-energy density) as an expansion in powers of the 
fields ip and a and the derivative operator d, with each 
term symmetric under the space group 25 (or a subgroup, 
if a perturbation explicitly breaks symmetries) . 

The terms involving only the matter fields ip are those 
that are gauge invariant, scalars under the space group, 
and time-reversal even. These can be found by consid- 
ering all combinations of the bilinears ^a^^ introduced 
in Section |V A[ and projecting into the identity repre- 
sentation. In the fully symmetric case, the only allowed 
term at first order in the bilinears is therefore given by 
the combination belonging to irrep Fi, labeled "+" in 
Table im 

The higher-order terms are found by taking symmet- 
ric combinations of all products of the bilinears. A gen- 
eral procedure exists to decompose a product into its 
irreps (see, for example, Ref. 44), but the projector for 
the identity representation can be constructed by sim- 
ply summing the transformation matrices over the space 
group. The set of all second-order combinations of the 
bilinears can be written as and the corresponding 

projector into the identity representation is given by an 
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expression analogous to Eq. (58), 



,(2) 



(61) 



lee 



where denotes an outer product. Subsequent terms 
in the expansion can be found by taking correspondingly 
higher-order products. 

For the fully symmetric case, the quartic terms in fact 
have a simple expression in terms of the bilinears as 



ri2 



(62) 



i.e., for each irrep F, a sum over the squares of the rows 
a. When perturbations reduce the symmetry, a similar 
result applies, but with the bilinears assigned to irreps 
of the appropriate subgroup of ©. This can be found by 



making use of a projector defined analogously to Eq. (61 1, 
but with the sum over the subgroup. 

The terms involving a and d are strongly constrained 
by gauge symmetry, Eq. ( 46 1 , and the lowest-order terms 



take the familiar form for a continuum gauge theory. 
They consist of a covariant-derivative term for each mat- 
ter field (with "minimal coupling" to a) and a pure gauge 
term, and can be viewed as deriving directly from the 
first two terms of the microscopic energy function Hj^ g 
in Eq. (10). The precise form of these terms can again 



be found by constructing the appropriate projection ma- 
trices, using the transformation properties of a and D 
specified in Section [IV E[ 

In the fully symmetric case, the long-wavelength action 
takes the form 

5 - ^(1 + 6%)(A<P«)nAv«) + X "I' 

ijX 

+ Rj2^l^jx + S''' + ---- (63) 

3X 

Terms involving more than two derivatives and matter 
fields beyond quartic order have been omitted, and are 
represented by the ellipsis. 

This model has a Coulomb phase for sufficiently large 
R, and can exhibit various types of Higgs transitions 
when R is reduced. The critical value Rc and the na- 
ture of the ordered phase will depend on the coefficients 
of the terms appearing in S'^^K Note that the action in- 
volves a total of 12 complex fields ipj^e and the constant 
quadratic term in fact has full SU(12) symmetry. While 
this is reduced considerably by S^'^\ permutation sym- 
metry remains between the 6 wavevectors n, and it is not 
possible in this case to eliminate any of the fields (pj ^e- 

Higher-order terms than those displayed in Eq. ( 63 ) 



are sometimes necessary to determine the ordering com- 
pletely. This occurs because the action including the 



quartic terms in Eq. ( 62 ) has continuous symmetry un- 



der rotations of the vector $ formed from the rows of 
a given irrep F (of dimension larger than one), while the 



microscopic model has only discrete spatial symmetries. 
There can therefore be states with different orientations 
of within the continuous manifold that are not re- 
lated by physical symmetries and that will be distin- 
guished by higher-order terms in the action. (A concrete 
example will be described below in the case with reduced 
symmetry.) In the case of a continuous transition, the 
higher-order terms are likely to be "dangerous irrelevant" 
(strictly irrelevant at the transition but important away 
from it) and the continuous symmetry is emergent at the 
critical point. 

Applying a perturbation coupling to one of the ob- 
servables adds to S" a first-order coupling to the corre- 
sponding bilinear $Jj . For those perturbations that break 
no symmetries, the appropriate coupling is to the bilin- 
ear for the identity irrep , and simply has the effect of 
changing the coefficient R. Symmetry-breaking pertur- 
bations, by contrast, lead to coefficients for the quadratic 
terms V*jxe'Pjxe thaX depend on j, x, and (for those that 
break translation symmetry) i. They can therefore dis- 
tinguish a subset of the fields (pjx£ as relevant at the 
transition, while the others can be eliminated from (in 
principle, integrated out of) the long-wavelength action. 
Besides reducing the symmetry-implied restrictions on 
terms in the action, applied perturbations can therefore 
reduce the number of matter fields that must be included. 

A particularly interesting case is where the symmetry 
is reduced to the extent that a single wavevector is rele- 
vant in the long- wavelength limit. According to Table |llj 
the bilinear corresponding to the observable (illus- 
trated in Figure [5]) is 

r+ 1 
^2'"' = (^^0,1+ - *o,i-) 



1 



($0,2- + $0,3- - $0,2+ - $0,3+) , (64) 



SO it is sufficient to add a coupling to di_ (with appro- 
priate sign) to pick out the bilinears with k = (It 
is of course possible for the same symmetries to be bro- 
ken spontaneously instead, in which case the observable 
di_ will become nonzero despite having no external cou- 
pling.) 

With such a perturbation, the unbroken symmetries 
Q are those for which k,q(k,i^) — and can be de- 

termined by inspection of the links in Figure [5j They 
include rotations by tt about all three cubic axes it^, and 
screw rotations by ^ along the helices formed by links 
with a given sign. There are no remaining improper 
symmetries, since these exchange chiralities and relate 
wavevectors Kjy^ with opposite x- 

The action is in this case given by 

Sr = + W,i)(A¥'i+)^(A¥'i+) + K\d X a|2 



(65) 



where other matter fields ^pj-^ have been eliminated 
(renormalizing coefficients such as i?), and coordinates 
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along Ml have been rescaled to make K isotropic, 
quartic terms can be written as 



The 



(66) 



where = = ip\^^cr^ip-^^^. Bihnears involving 

matter fields other than (pi_|_ have been dropped, and the 
equality between ui and U2 previously enforced by rota- 
tion symmetry no longer applies. Because of the reduced 
symmetry, certain observables such as spin spirals other 
than £1+ can no longer be expressed in terms of relevant 
matter fields; these cannot therefore be order parameters 
for any transitions described by the action S'r. Referring 
to Table [nj one also finds that the observables — and 
couple to $0 = '^'o,i+,i+ and so are always nonzero 
(breaking no further symmetries). 

The action 5^ is similar to that describing transitions 
of the cubic dimer model,'2iH23l having two relevant mat- 
ter fields and quadratic terms that are symmetric un- 
der SU(2). Several symmetry-breaking scenarios in this 
model have been considered by Chen et al.,1^ and use 
will be made of their results at several points in the 
following. (The same critical theory also describes the 
Mott insulator-superfluid transition for quantum bosons 
at half filling on a square lattice.—) 



B. Ordered states 



In the models of Eqs. (63) and (65), a transition out of 



the Coulomb phase occurs when the matter field con- 
denses and the gauge field <y. acquires a finite correlation 
length by the Higgs mechanism. In both cases (f has 
multiple components, so condensation also breaks cer- 
tain spatial symmetries. Because of the gauge symmetry 
of if, it is more convenient to describe these transitions 
in terms of the gauge-invariant bilinears which will 
therefore serve as order parameters. (It is possible to ap- 
ply sufficient perturbations that there is a single relevant 
field ^p remaining, so there is no spontaneous symmetry 
breaking. An example will be treated briefly at the end 
of this section.) 

In the fully symmetric case, the Higgs transition in- 
volves the matter fields condensing at some subset of the 
6 wavevectors Kj^, with an orientation chosen for each 
two-component complex vector y>^-^ (up to gauge redun- 
dancy). This clearly admits the possibility of a direct 
transition where the matter field condenses into such a 
configuration. One could instead imagine a transition 
at which the Coulomb-phase correlations remain, but 
spatial symmetry is spontaneously broken, followed at a 
lower temperature by a Higgs transition of the remaining 
relevant matter fields. The latter will be identical to one 
with explicitly broken spatial symmetry, to be discussed 
below, and may be continuous or of first order. 

An interaction between second-neighbor spins, as in 
Figure [4] provides a specific example of a perturbation 
that preserves all symmetries but can lead to an ordered 




FIG. 7: Second-neighbor interactions between spins, as in 
Figure |4] favoring a spin spiral, as in Figure [5] Such a spiral, 
one of 12 symmetry-related states, minimizes the energy if 
all second-neighbors have antiferromagnetic (AF) Heisenberg 
interactions (so no symmetries are explicitly broken). The 8 
solid (black) lines show second neighbors that satisfy this per- 
turbation, while the 4 shown dashed (and red) are unsatisfied. 
Considering all configurations of two adjacent tetrahedra (of 
6 pairs, at most 4 are AF aligned) shows that one cannot sat- 
isfy more than 2/3 of the second- neighbor links while obeying 
the ice rule. (It is worth noting a counterintuitive result of 
the geometry of spin ice: second neighbors that are both AF 
aligned with their common neighbor point towards tetrahe- 
dra of the same orientation and so are also AF aligned with 
each other.) Conversely, ferromagnetic interactions between 
second-neighbor spins disfavor the spiral configuration, and 
can instead lead to a bond-ordered state such as shown in 
Figure [8] 



state. It is demonstrated in Figure [7] that antiferromag- 
netic interactions are minimized by arranging the spins 
into a spiral state (as in Figure [5|, which is one of 12 
symmetry-related configurations. (The argument does 
not prove that these 12 are the only minimal-energy con- 
figurations; additional interactions may in practice be 
necessary to stabilize the spiral state.) A spiral with 
axis along Uj and chirality x involves condensation of 
the matter field at wavevector Hj^] the diagonal bilinear 



3,3X,3X 



3X 



corresponding to the observable 
£jp^, serves as an order parameter and can be positive or 
negative in the ordered phase. 

In cases where a perturbation reduces the symmetry 
explicitly, the types of ordered states that are preferred 
will again depend on the details of the perturbation. As 
long as there is sufficient symmetry that more than one 
wavevector k, is relevant to the ordering, one has a similar 
situation to the fully symmetric case. 

The case where the single wavevector is selected 
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has been considered above: the action is given in Eq. ( 65 ) 
and a suitable perturbation is iUustrated in Figure [5J 
Also shown is a configuration of spins that minimizes 
the energy and may be selected by an ordering transi- 
tion. The spins form a right-handed spiral described by 
the observable £1+ ^ $1 = ip\^cT^ip^j_, using the nota- 
tion of Eq. (66). There arc two such states, related by 
time reversal (or the translation T2, for example), with 
opposite sign for <i>i. 

One can therefore expect a transition in the presence of 
this perturbation into a Higgs phase with a nonzero value 
of the real order parameter $1. (One of the two linear 
combinations ^i+fl i (ySi+.i condenses.) This transition, 
at which an Ising symmetry is broken, resembles that in 
the "2-GS" model of Chen et al.,!^ which was found to 
be of first order. 



The form of the quartic terms in Eq. ( 66 1 suggests the 



possibility of other ordering transitions, where the vector 
# aligns elsewhere in the 3-dimensional space. Accord- 
ing to Eq. ( 56 1 , the components $2 and $3 are even un- 
der time reversal, so states where only these are nonzero 
are paramagnetic; the appropriate order parameters can 
be found by reference to Table |llj Changing the relative 
magnitude of ui and U2 tunes the system between these 
and the spin-ordered states with ($1) 7^ 0, and can be 
effected by further perturbations that do not reduce the 
symmetry but disfavor spin-ordering. An example is pro- 
vided by the second-neighbor interactions illustrated in 
Figure |4] the argument of Figure [7] implies that the spin 
spiral maximizes the energy when these are ferromag- 
netic. 

An example of a state with (<i>2) ^ and ($1) — 
(<i>3^ is illustrated in Figure |8] According to Ta- 
ble O (in particular, the rows for irreps and X^), 
$2 = ^2,1-1- can be associated with the observables 
'^2o(~f'"2) and d2o(~f'"2): which transform identically 
under the reduced symmetry group. In the ordered state 
shown, links with additional antiferromagnetic interac- 
tions are fully satisfied only in alternating (010) planes, 
so the order parameter has wavevector — f ■U2- For tetra- 
hedra with cr = (orange), 6,2 is nonzero and staggered 
on successive (010) planes, so $2 '^2o(— §^2) 7^ 0. By 
contrast, ^2 vanishes on tetrahedra with a = 1 (blue), 
so $3 ~ d2i(— §1*2) — 0. (The single spin configuration 
shown is of course only a member of the ensemble; the 
state has no spin order.) 

This phase resembles the intermediate paramagnets 
found by Chen et ahl^in the cubic dimer model. The lat- 
ter occur when sufficient perturbations are applied that a 
Higgs transition occurs without spontaneously breaking 
any symmetries, followed by a lower-temperature tran- 
sition into an ordered phase. The bond-ordered phase 
illustrated in Figure [8] does break spatial symmetries at 
the Higgs transition but preserves time-reversal symme- 
try and has no spin order. (It is also likely an intermedi- 
ate phase: reducing the temperature further will remove 
the residual entropy leading to a spin-ordered state such 
as in Figure [5]) 




FIG. 8: Illustration of the bond-ordered phase, with o rder p a- 



VI B 



to 



rameter $2,i+,i+ ~ oi2o(— f W2), proposed in Section 
occur in the presence of the perturbation shown in Figure [5] 
Only those links with additional antiferromagnetic (AF) in- 
teractions are shown (in red), and in this putative ordered 
state these links are divided into two types by their spin cor- 
relations. The thicker solid lines show links that are strongly 
AF aligned, i.e., pairs of spins with large negative expecta- 
tion value for their inner product, while the thin dashed lines 
indicate weaker AF alignment. (If all the antiferromagnetic 
bonds were fully satisfied, the state would freeze into the spin 
spiral shown in Figure [s]) Symmetry under the translation 
T2 and the remaining ^ screw rotations (along ui) are bro- 
ken, making the two tetrahedron orientations (labeled "0" and 
"1") inequivalent. The spins in the picture show an extreme 
example of such a configuration, where all favored antiferro- 
magnetic bonds are satisfied and the remainder are unsatis- 
fied; the actual ordered state has no spin ordering and a finite 
spin correlation length in all directions, including along the 
favored chains. 



The transition, by contrast, is apparently in the same 
class as the "4-GS" model of Chen et al.,"^- described by a 
matter-field doublet and easy-plane anisotropy for the 
bilinear Unlike the case where $1 becomes nonzero 
and a single component of the complex vector ^Pi^ con- 
denses, there is now a connected low-energy manifold, 
described by the phases of the two components of Vi^. 
This led Chen et al.l^ to conjecture that a continuous 
transition is possible, although a weakly first-order tran- 
sition was found in their model. 

Note that this example illustrates the observation 
made in Section [VI A| that higher-order terms are some- 
times needed to specify the ordering. The quartic terms 
have a symmetry under continuous rotations mixing the 
$2 and 4>3 components of the vector (The quadratic 
term ip\^ipi^ = $0 = |*| is clearly also symmetric.) For 
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U2 < ui, the free energy is minimized for $ lying any- 
where in this plane, but this apparent degeneracy will 
be broken by higher-order terms such as <1'2'^3, allowed 
by the discrete microscopic symmetries. (These higher- 
order terms also prevent the appearance of unphysical 
Goldstone modes due to the broken rotation symmetryl^Sl 
All correlation functions are short-ranged in the ordered 
phase because the physical symmetries are discrete.) 

A particularly interesting special case occurs when the 
quartic coefficients ui and 1*2 are tuned to the same value. 
On the Higgs side of the transition, the order changes 
from one type to the other at this point, presumably in 
a strongly first-order "spin flop". At the multicritical 
point where this transition line meets the boundary of 
the Coulomb phase, the long-wavelength action given in 
Eq. (65 1 has an emergent SU(2) symmetry. There is then 
the possibility of a continuous transition in the "noncom- 
pact CP^" class believed to describe ordering in the cubic 
dimer model .l^iHlI] 

A final example, where theoretical arguments strongly 
suggest a continuous transition, is when an applied per- 
turbation favors a unique configuration such as the spin 
spiral shown in Figure [5j A suitable choice is clearly 
given by an applied magnetic field, coupling to the spins 
via a Zeeman term, with the same spatial configura- 
tion as the spiral. This perturbation breaks transla- 
tion symmetry under T2 and T3 (as well as time rever- 
sal), so not only picks out a single wavevector, but also 
breaks the degeneracy in £. The critical theory can there- 
fore be written in terms of a single scalar matter field 
ip = ^((^i+,o + ^-iid takes the standard form for 

a U(l) gauge theory. The Higgs transition breaks no spa- 
tial symmetries (($1) > 0, even in the Coulomb phase), 
and is in the usual inverted-XY universality class:^^ 



VII. CONCLUSIONS 

This work has analysed phase transitions out of the 
Coulomb phase in spin ice, and shown that they can be 
understood as Higgs transitions of an underlying gauge 
theory. A long-wavelength description results from find- 
ing the relevant modes of an emergent matter field, which 
transform according to representations of the magnetic 
symmetry group. The possible Higgs transitions are de- 
termined by forming order parameters from the matter 
fields and studying an effective action written in terms 
of these and the fluctuating gauge field. This stands in 
contrast to the standard paradigm of phase transitions, 
where the Ginzburg-Landau action is expressed directly 
in terms of the order parameters.^^fil 

The phenomenology is particularly rich in the case of 
spin ice, as a consequence of symmetry under inversion 
of all spins ("time reversal"), which, as argued in Sec- 
tion 
two 



transition involving these wavevectors seems quite un- 
likely. The long-wavelength theory is instead expressed 
in terms of matter fields at 6 wavevectors, leading to 
a variety of possible ordering transitions, both in the 
fully symmetric case and when applied perturbations re- 
duce the number of relevant modes of the matter field. 
Because of this complexity, no attempt has been made 
here to describe all possible transitions or ordered states, 
and concrete results have been provided only in cases 
of particular interest. The analysis presented here can 
nonetheless be adapted to study other transitions in this 
system that may be discovered by experiments, numerics, 
or other theoretical approaches. 

In fact, it is interesting to note that the Monte Carlo 
studies of Melko et al.l^ found that the ground state in 
the presence of dipolar interactions is given by exactly 
the type of spin spiral shown in Figure [s] (As argued in 
Figure [7] this spiral is also favored by antiferromagnetic 
Heisenberg interactions between second neighbors.) The 
first-order transition observed in the numerical simula- 
tion is interpreted in the language of the present work as 
a Higgs transition of the matter field ip^^ such that the 
diagonal bilinear ^jjxJx acquires a nonzero expectation 
value. 

One straightforward prediction of the present work is 
a continuous transition in the inverted-XY universality 
clas^^ into the spin spiral shown in Figure [sj in the 
presence of a helical magnetic field with the same spatial 
structure.!^ While this transition does not spontaneously 
break any symmetries, it is nonetheless an example of an 
exotic "non-Landau" transition, described by a critical 
theory of a single matter field coupled to an emergent 
U(l) gauge field. Producing a static magnetic field with 
the required structure on the atomic scale would of course 
be infeasible in an experiment, but trivial in a numerical 
simulation. 

An intriguing possibility that arises naturally from the 
analysis of Section VI B is an emergent SU(2) symme- 
try at the point shown in Figure [9] where the Coulomb 
phase meets two distinct ordered phases. The transitions 
into the two low-temperature phases, with spin-spiral and 
bond-based order parameters respectively, resemble those 
studied by Chen et al.l^in the cubic dimer model, which 
were shown to be of first order. At the point where the 
three phases meet, however, the long-wavelength theory 
given in Eq. (65) has an emergent SU(2) symmetry, al- 



lowing a continuous transition described by the noncom- 
pact CP^ (NCCP^) critical theory. If such a multicriti- 
cal point is indeed found in numerics, its properties can 
be compared to those for other transitions^^ '^'^ supposed 
to be described by NCCP^. A clear signature of such 



IV B leads to an extra (Kramers) degeneracy at the 



lighest-symmetry wavevectors. Since they cannot be 
dispersion minima of any hopping Hamiltonian (preserv- 
ing time-reversal symmetry; see Appendix^, a Higgs 



a scenario would be an emergent SO (3) symmetry^ in 
the correlations of the components of the vector order 
parameter as described in Section [Vl B[ 

It must be noted that, even in these cases, the pos- 
sibility of a first-order transition cannot be excluded on 
the basis of the analytic arguments presented here. Al- 
though the long-wavelength theory may exhibit a con- 
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FIG. 9: Schematic phase diagram as a function of tempera- 
ture T and antiferromagnetic interaction V2 between second- 
neighbor spins, in the presence of a (fixed) perturbation such 
as shown in Figure [5] The spin spiral, also shown in Figure [5] 
and the bond-ordered phase. Figure |8j are ordered phases 
with distinct order parameters and are separated by a strongly 
first-order ("spin fiop") transition. The transitions into the 
Coulomb phase from the spin spiral and bond-ordered phase 
are described in Section fVI B[ At the multicritical point where 
the three phases meet, there is an emergent SU(2) symmetry, 
allowing a continuous transition described by the noncompact 
CP^ theory. 



tinuous transition (which remains controversiaPSESl for 
NCCP^), this can be preempted by a first-order transi- 
tion in a given microscopic model. In such cases, it may 
nonetheless be possible to tune the latent heat to zero us- 
ing additional perturbations, as has been demonstrated 
explicitly in the cubic dimer modebl^^^ 

The analysis presented here excludes states with 
nonzero net polarization, which are favored whenever a 
uniform magnetic field is applied.l^^^i^ (The "squiggle" 
state found to occur in quantum ic^^ also has nonzero, 
though unsaturated, polarization. Quantum effects can 
be included within the present approach by the addition 
of interactions favoring "flippable" configurations.) To 
remove this restriction, one must include the possibility 
of nonzero uniform flux in Eq. (|6|, which has the effect 
of modifying the background gauge field a. Results for 
this case will be presented elsewhere. 

This work has also treated only the case where the ice 
rule is uniformly obeyed, and Eq. ^ can be viewed as 
a constraint. This is equivalent to forbidding magnetic 



monopoles,'^ and requires temperature much lower than 
the scale of the nearest-neighbor interactions. Deviations 
from this limit can be treated at a simple level by anal- 
ogy with finite-size scaling: the analysis presented here is 
strictly applicable only in the thermodynamic limit with 
zero monopole density, but has relevance even when the 
system size and mean monopole separation are finite but 
large. 
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Appendix A: Calculation of MSG commutation 
relations 



The commutation relation, Eq. (21 1, of the FCC trans- 



lation operators 7i and T2 was calculated in SectionpIIB| 



by application of the relation Eq. ( 18 ). Commutation re- 
lations for other operators Q can be found by the same 
approach, with the subtlety that these depend on the 
arbitrary choices of global phase made for the functions 

TTQ. 

For example, one has MiTi ~ T2M1 (where Mi is a 
proper rotation by 27r/3 about the diamond link 5i, with 
corresponding MSG operator TZi), so the relative phase 
of Ti and T2 can be chosen such that TZiTi = T2R-i- In 
other words, one fixes 



TT-Ri (TijVri (j) = TTr^ (Ml j)7rKi [j) 



(Al) 



by appropriate choice of the relative phase of ttj--^ and 
ttt-^. Using Eq. (18 1, this gives 



(j>f,(Mij)=c.*(j^Tij)c.(Mij 



iTTij), (A2) 



which determines other commutation relations involving 
71 and 72. One can similarly fix the phase of 73 relative 
to 7i and 72 by requiring TZiT^ = TiTZi- 

While these relations fix the relative phase of the trans- 
lation operators %, their absolute phases can be fixed 
using commutation relations such as T1M2 — M2T^^Ti. 
The corresponding expression in terms of the MSG oper- 
ators depends on the absolute phase of 72, as confirmed 
by considering 



^n-^rr^n^Tf^T, ij) = ^K.(J>f,(IR2J>7^.(T2"'TlJ>f^(T2-lTu>^^(J) 

= ^f^(T^iTij)c^(Tij ^ T1M2J ^ M2J ^ J ^ T^'Tij) , 



(A3) 
(A4) 



where Eq. ( 18 ) has been used to eliminate n-ji^ and ttj--^ . Sinc e TZ^ ^7?,2 7^ ^71 = ±£, this expression cannot depend 
on J, and so one can simplify the chain appearing in Eq. (A4| by making a suitable choice of j. It is particularly 
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convenient to set j = Tj^ "'^TjO, where 0, the origin of the coordinate system, obeys M2O = 0. In this case, one has 
simply 

^n^^rr^n.Tf^r, (Tr'T20) - 7rf^(0)^(T20 ^ 0) . (A5) 

Since the absolute global phase of ttj-^ is as yet unspecified, one can choose Trj-2{0)uj{0 — > T2O) = 1, giving 
'^n~^T^^n T^^T ~ ^' '^^ equivalently 71 72.2 = 7?.27^~^7i. 
Note that the quantity 

<^r.O') = 7r7-.(i)c^0'->T,j) (A6) 



is invariant under changes of the gauge of the background field a, with ttq transforming to maintain Eq. (16). The 
choice made after Eq. ( |A5[ ) gives 0X2(0) = 1; the same is true for Ti, according to Eq. (A2|, and for Ta similarly. 
Furthermore, Eq. ( 18 1 implies 

U) = <P% (O)c^(O ^ • • • ^ J ^ T, j ^ • • • ^ T,:0 ^ 0) ; (A7) 

the chain appearing on the right-hand side is a closed loop ('parallelogram') consisting of two pairs of parallel segments. 
Counting the number of enclosed hexagons gives 



(^■) = (_l)i("-'-.-.), (A8) 



where rj is the position vector of site j and aj G {0, 1} is its sublattice. Note that this quantity is real for all sites j 
of the diamond lattice. 

The commutation relation of a primitive translation % and any other MSG operator Q can be found by first 
considering the corresponding expression for the space-group operators and Q. It is simple to show that the 
combination QT^Q"^ is a pure translation T-^^^f,., by displacement MQe^. One is then led to consider the combination 

7i^;Je,2'^2~^ and to compute 

^r„^ Qr.Q-(j)-4- (j)7r^,(Q-i)c.*(T,Q-ij^Q-^j>(TM^„e,j^j), (A9) 



where Eq. (18) has been used to eliminate ttq. Both 



chains in this expression have total displacement equal to 
a minimal FCC vector and can consistently be taken to 
follow the conventional path, as in Eq. (20). One there- 
fore finds 



(j) = '^fM„..(j>r.(Q-'j) 



(AlO) 



which is valid for all sites j, but, since Tj^^^^ QT^Q ^ = 
E, must in fact be independent of j. Taking j = gives 
''Q-iQ = "Mq"'^ Vq, and one can choose the sign of Tmqei 
to make (?!>rM,jje, (0) = 1, so 

QT^ = (-l)5('^«-^«"'+"^VMQe.Q, (All) 



using Eq. ( A8 ) 



The MSG cubic translation K-i defined in Eq. ( 22 ) can 
alternatively be expressed as 



(A12) 



where [i ± 1] appearing in the subscript indicates arith- 
metic modulo 3. Analogous to 07-, defined in Eq. (A6I, 
the quantity 

(A13) 



is gauge invariant. The chain uj{j —i' Kij), consist- 

ing of minimal FCC translations interpreted according 
to the convention used in Eq. ( |20| ) , is chosen to comprise 
a right-handed helix as in Figure [3] Expressing tt/c^ in 
terms of nj-. with Eqs. ( 19 ) and ( |A12[ ) and then using 
Eq. ([A7|, this gives 



0/c. U) 



-c^T,.^,, (T[,+i]T-V)0f.(T[,+i]Tr^j>7].+,j (j) 

(A14) 

1 , (A15) 



the result expressed as Eq. ( 23 1 in Section III B 1[ 



Appendix B: Microscopic hopping model 

A significant part of the present work involves finding 
representations of the magnetic symmetry group (MSG), 
to which the important modes of the matter field '0 be- 
long. Certain results of this analysis can be confirmed 
by studying a solvable model with the same spatial sym- 
metries and finding the degeneracies and transformation 
properties of its modes. An obvious choice is to use the 
analogy, noted in Section [llj with a quantum model of 
noninteracting bosons on the diamond lattice in the pres- 
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ence of a static gauge field a\ . In this appendix, tfie spec- 
trum of this (three-dimensional) Hofstadter problem'^^ 
will be presented. 

Starting from the model of Eq. (10), softening the 
unit-modulus constraint on tp — e'^ and dropping the 
fluctuating gauge field A\ leaves a model with the ap- 
propriate symmetries and hopping terms only between 
nearest-neighbor sites. It is nonetheless worthwhile to 
include further-neighbor hopping, and to study a generic 
model with the same symmetries. As will be noted be- 
low, the nearest-neighbor problem has additional "acci- 
dental" degeneracies that cannot be expected to survive 
in the presence of the fluctuating gauge field. 



1. Hopping phases 



Time reversal Q acts according to Eq. (27), and is a sym 



metry of the hopping Hamiltonian only if ttq can be found 
so that 



(B4) 



for all sites i and j. Taking the product of both sides 
for the links i ^ j around a hexagonal plaquette tt, 
this implies that J^(j_j.j>)g^ tij must be real if Q is un- 
broken. This is indeed the case for the phases a;(i — 
j) = e^'^'^'^j , for which the product is always —1. Per- 
turbations that break Q (such as an applied magnetic 
field) are equivalent in the hopping problem to a non- 
trivial magnetic flux (not zero or half a flux quantum) 
through the hexagonal plaquette. 



The hopping model is written in terms of complex 
fields tpi and has Hamiltonian 



(Bl) 



where the (dual diamond) sites i and j need not be neigh- 
bors and every pair is counted twice. For any space-group 
transformation Q, the hopping amplitudes Uj transform 
under the corresponding MSG operator Q, and so obey 



> j)uj*( 



— ^ 



(B2) 



where Eq. ( 18 ) has been used. In order for the energies 
to be real, one requires — t* 



For nearest-neighbor links, hopping phases obeying 
Eq. (|B2t can be extracted directly from the first term 



of Eq. (10): for these links, take ti 



t, 



\tij\uj*{i — ?> j). A consistent set of phases for the further- 
neighbor hopping is given by 



(0) 



tij — 



(B3) 



where t^^^ is real and is a complete set of symmetry- 
related paths from i to j. In other words, start by choos- 
ing an arbitrary path Ci = i — > • • • ^ j, and sum over all 
space-group operators Q such that QCi also runs from 
i to j (i.e., Qi = i and Qj — j). It is then straightfor- 



magnitudes i,-^'* are symmetric under the space group. 
For certain pairs of sites, the sum in Eq. (B3) van- 



ward to show that Eq. (B2 ) is satisfied, provided that the 



ishes and the hopping tij must vanish. This occurs, for 
instance, when the two sites are at opposite ends of a 
hexagon of the diamond lattice. The set €ij can be 
taken as consisting of the two shortest paths around the 
hexagon, but, as noted in Section pIIB[ these have u dif- 
fering by a minus sign, so the sum vanishes. 

To treat the case where perturbations reduce the sym- 
metry of the problem, one can modify the hopping am- 
plitudes to obey Eq. (B2) for the appropriate subgroup. 



2. Spectrum 

Once the hopping phases are determined, the Hamil- 
tonian TJhopping can be expressed in the basis of eigen- 
states of ICi. The structure of the energy eigenstates 
within the cubic unit cell can be found from the effective 
K-dependent hopping 



t. 



,-(/«) E^Kij7rK;(i)e 



iRK 



(B5) 



where i and j are both sites within the cubic unit cell of 
the diamond lattice. The sum is over all cubic transla- 
tions K by displacement iijj- 

The dispersion for the fully-symmetric case, shown in 
Figu re \T0\ exhibits various features anticipated in Sec- 
tion |IV[ Every mode is twofold degenerate, labeled by 
£ = 0,1, and the energies at wavevectors k and k,q{k,) are 
equal for all Q. In particular, inversion symmetry implies 
degeneracy between the states at k and [k,— jSq], such as 
F and R. As noted in Section [IV B[ there are Dirac cones 
in the spectrum at the latter two points. The global min- 
ima of the spectrum are at the points X and M, which, as 
noted in Section IV C[ form a star containing 6 wavevec- 
tors. (With only nearest-neighbor hopping, the energy is 
equal along the line Z between these points, but this acci- 
dental degeneracy is lifted by further- neighbor hopping.) 

The symmetry-breaking perturbation introduced in 
Section VI A allows a single wavevector Ki+ to be distin- 
guished from the rest. This can be confirmed by calculat- 
ing the spectrum with an additional hopping term that 
breaks the same symmetries as the perturbation shown 
in Figure [5] With the appropriate choice of sign for this 
hopping, there is indeed a unique minimum energy at the 
wavevector 

A perturbation that selects a single spiral state must 
break time-reversal symmetry as well as the FCC transla- 
tions Ti, and corresponds to additional hopping phases in 
this model. The applied magnetic field favoring the spin 
order shown in Figure [5] can be produced with a vector 
potential having the same spiral structure and the same 
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cubic unit cell as the ordered state itself. The resulting 
spectrum is similar to the previous case, with a unique 
minimum at but with the twofold degeneracy also 
broken. 



FIG. 10: Mode energy (in units of the nearest-neighbor hop- 
ping amplitude) in a hopping model with the full symme- 
try of the diamond lattice, as a function of wavevector k 
within the cubic reduced Brillouin zone OSr. There are 8 
sites within the cubic unit cell, and hence 8 bands, but there 
is a twofold degeneracy due to the anticommutation of the 
FCC translations Ti. As discussed in Section [IV B[ there are 
additional degeneracies at the points V and R, which are each 
surrounded by (three-dimensional) Dirac cones, visible here 
as crossings. The dashed line shows the dispersion with hop- 
ping only between nearest-neighbor sites, which has fiat bands 
along the lines Z joining points X and M. (Throughout this 
appendix, all labels for symmetry points are referred to 58r.) 
For nonzero dispersion along this line, one requires hopping 
between at least sixth neighbors; the solid line has a weak 
additional hopping between sites separated by 2ei +62 — 63 
and equivalent vectors. 



M. J. Harris, S. T. Bramweh, D. F. McMorrow, T. Zeiske, 
and K. W. Godfrey, Phys. Rev. Lett. 79, 2554 (1997). 
S. T. Bramwell and M. J. P. Gingras, Science 294, 1495 
(2001). 

P. W. Anderson, Phys. Rev. 102, 1008 (1956). 

L. Pauhng, J. Am. Chem. Soc. 57, 2680 (1935). 

A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, 

and B. S. Shastry, Nature 399, 333 (1999). 

T. Fenneh, O. A. Petrenko, B. Fak, S. T. Bramweh, M. 

Enjalran, T. Yavors'kii, M. J. P. Gingras, R. G. Melko, 

and G. Balakrishnan, Phys. Rev. B 70, 134408 (2004). 

R. G. Melko and M. J. P. Gingras, J. Phys.: Condens. 

Matter 16, R1277 (2004). 

T. Fenneh, P. P. Been, A. R. Wildes, K. Schmalzl, D. Prab- 
hakaran, A. T. Boothroyd, R. J. Aldus, D. F. McMorrow, 
and S. T. Bramweh, Science 326, 415 (2009). 
C. L. Henley, Annual Review of Condensed Matter Physics 
1, 179 (2010). 

C. Castelnovo, R. Moessner, and S. L. Sondhi, Nature 451, 
42 (2008). 

R. Youngblood, J. D. Axe, and B. M. McCoy, Phys. Rev. 
B 21, 5212 (1980). 

D. A. Huse, W. Krauth, R. Moessner, and S. L. Sondhi, 
Phys. Rev. Lett. 91, 167004 (2003). 



C. L. Henley, Phys. Rev. B 71, 014424 (2005). 

S. V. Isakov, K. Gregor, R. Moessner, and S. L. Sondhi, 

Phys. Rev. Lett. 93, 167204 (2004). 

S. V. Isakov, R. Moessner, and S. L. Sondhi, Phys. Rev. 
Lett. 95, 217201 (2005). 

R. Moessner and S. L. Sondhi, Phys. Rev. B 68, 064411 
(2003). 

^"^ S. Yoshida, K. Nemoto, and K. Wada, J. Phys. Soc. Jpn. 
73, 1619 (2004). 

L. D. C. Jaubert, J. T. Chalker, P. C. W. Holdsworth, and 
R. Moessner, Phys. Rev. Lett. 100, 067207 (2008). 
^'^ S. Poweh and J. T. Chalker, Phys. Rev. B 78, 024422 
(2008). 

2° L. D. C. Jaubert, J. T. Chalker, P. C. W. Holdsworth, and 
R. Moessner, Phys. Rev. Lett. 105, 087201 (2010). 
A. Andreanov, J. T. Chalker, T. E. Saunders, and D. Sher- 
rington, Phys. Rev. B 81, 014406 (2010). 
R. G. Melko, B. C. den Hertog, and M. J. P. Gingras, Phys. 
Rev. Lett. 87, 067203 (2001). 

2^ T. S. Pickles, T. E. Saunders, and J. T. Chalker, Europhys. 
Lett. 84, 36002 (2008). 
C. Xu, Phys. Rev. B 81, 144430 (2010). 
P. W. Anderson, Phys. Rev. 130, 439 (1963). 

^® L. D. Landau and E. M. Lifshitz, Statistical Physics, 



23 



Butterworth-Heinemann, New York (1999). 
'^'^ O. I. Motrunich and A. Vishwanath, Phys. Rev. B 70, 
075104 (2004). 

T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and 
M. P. A. Fisher, Science 303, 1490 (2004); T. Senthil, L. 
Balents, S. Sachdev, A. Vishwanath, and M. P. A. Fisher, 
Phys. Rev. B 70, 144407 (2004). 

L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and K. 

Sengupta, Phys. Rev. B 71, 144508 (2005); Prog. Theor. 

Phys. Suppl. 160, 314 (2005). 
^° F. Alet, G. Misguich, V. Pasquier, R. Moessner, and J. L. 

Jacobsen, Phys. Rev. Lett. 97, 030403 (2006). 
^1 S. Powell and J. T. Chalker, Phys. Rev. Lett. 101, 155702 

(2008); Phys. Rev. B 80, 134413 (2009). 

D. Charrier, F. Alet, and P. Pujol, Phys. Rev. Lett. 101, 

167205 (2008). 

G. Chen, J. Gukelberger, S. Trebst, F. Alet, and L. Ba- 
lents, Phys. Rev. B 80, 045112 (2009). 
'^^ S. Papanikolaou and J. J. Betouras, Phys. Rev. Lett. 104, 
045701 (2010). 

D. Charrier and F. Alet, Phys. Rev. B 82, 014429 (2010). 
D. L. Bergman, G. A. Fiete, and L. Balents, Phys. Rev. B 
73, 134402 (2006). 

S. Sachdev, Quantum Phase Transitions, Cambridge Uni- 
versity Press, Cambridge (2011). 
^® N. Shannon, O. Sikora, F. PoUmann, K. Penc, and P. 
Fulde, arXiv:1105.4196 (unpublished). 
Z. Nussinov, C. D. Batista, B. Normand, and S. A. Trug- 
man, Phys. Rev. B 75, 094411 (2007). 
M. E. Fisher and A. Aharony, Phys. Rev. Lett. 30, 559 
(1973). 

D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 

''^ One can instead justify the introduction of the charges 
using a dual description in terms of an effective electric 
field. Because the magnetic field is discrete, the electric 
field is compact (defined only modulo 2tt) and hence al- 
lows monopoles. These electric monopoles are exactly the 
matter fields described by 6. 

J. F. Cornwell, Group Theory in Physics: An Introduction, 
Academic Press (1997); Group Theory in Physics, Vol. I, 
Academic Press (1984). 

M. Lax and J. J. Hopfield, Phys. Rev. 124, 115 (1961). 

E. Brown, Bull. Am. Phys. Soc. 8, 256 (1963); Phys. Rev. 
133, A1038 (1964). 



J. Zak, Phys. Rev. 134, A1602; 134, A1607 (1964). 
■'■^ H. Overhof and U. Rossler, Phys. Stat. Sol. 26, 461 (1968). 
This is sometimes known as the "projective symmetry 
group" since it provides a projective representation of 
the space group. 

X.-G. Wen, Phys. Rev. B 65, 165113 (2002). 
The multiplication rules for the MSG operators are equal 
to those of the space-group operators up to a sign, so the 
group consisting of ±Q is by definition a "double group" ^ 
(Such groups appear in the context of fermions, where cer- 
tain rotations change the sign of the spinor wavef unction.) 
The relevant representations are those where the multipli- 
cation rules have the same sign as those of the MSG oper- 
ators, the so-called "extra" representations of the double 
group. The analogy with the fermionic case is not com- 
plete, however, since the multiplication table is different 
here from the fermion case (where, for example, all trans- 
lations commute), and so tabulated extra representations 
are of no use. 

^'^ C. Dasgupta and B. L Halperin, Phys. Rev. Lett. 47, 1556 
(1981). 

This transition can also be understood using the classical- 
to-quantum mapping of Refs. [18| and 1191 The nonuniform 
magnetic field maps to a potential energy distinguishing 
the two sublattices of the resulting quantum model, and 
the spin spiral becomes a Mott insulator. The quantum 
phase transition from the superfluid (corresponding to the 
Coulomb phase) is in the standard 2 -I- ID XY univer- 
sality class, in agreement with the present conclusions. 
(The transition is "inverted" because the superfluid, the 
ordered phase of the quantum model, maps to the high- 
temperature Coulomb phase of the classical model.) 
This resembles the SO (5) scenario proposed for the 
cuprates,!^ where two apparently unrelated order param- 
eters form a vector with an approximate symmetry. In the 
present case, there is an exact emergent SU(2) symmetry 
precisely at the multicritical point, provided that higher- 
order terms are (as seems likely) irrelevant. 
^'^ E. Demler, W. Hanke, S.-C. Zhang, Rev. Mod. Phys. 76, 
909 (2004). 

A. B. Kuklov, M. Matsumoto, N. V. Prokof'ev, B. V. 
Svistunov, and M. Troyer, Phys. Rev. Lett. 101, 050405 
(2008). 



